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Evaluation  of  the  original  three -component  adaptive  processor 
revealed  that  improvements  in  signal -to -noise  ratio  are  much  greater  for 
Rayleigh  waves  than  for  Love  waves.  A new  algorithm  has  been  proposed  for 
Love  wave  processing  to  put  the  effectiveness  of  the  Love  wave  processor  on 
a par  with  that  of  the  Rayleigh  wave  processor.  This  new  algorithm  tracks 
the  incoming  Love  in  azimuth  and  passes  its  entire  duration.  In  addition, 
energy  from  azimuths  outside  some  pre-set  limit  is  rejected  in  a form  of 
azimuthal  filtering. 

In  this  report,  expected  signal-to-noise  ratio  gains  over  the 
bandpass  filter  are  determined  for  both  the  original  and  new  versions  of  the 
three-component  adaptive  processor.  Detection  threshold  improvements  are 
then  determined  by  applying  both  versions  to  two  event  suites.  The  evaluation 
is  performed  for  both  single-site  and  beam  data. 
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SECTION  1 
INTRODUCTION 


The  three-component  adaptive  processor  (TCA)  was  developed 
at  the  Lamont  Geophysical  Observatory  (Shimshoni  and  Smith,  1964)  and  was 
evaluated  for  single-site  and  beam  data  by  Texas  Instruments,  Inc.  (Lane, 
1973).  This  processor,  designed  to  improve  the  detectability  of  long-period 
Rayleigh  and  Love  waves,  takes  advantage  of  the  known  phase  relationships 
among  the  three  mutually  perpendicular  long-period  seismometer  traces. 
Improvements  in  the  signal-to-noise  ratio  (SNR)  can  be  achieved  when  these 
phase  relationships  are  utilized  in  the  processor  design. 

The  evaluation  of  the  original  TCA  processor  (Lane,  1973) 
revealed  that  improvements  in  signal-to-noise  ratio  are  much  greater  for 
Rayleigh  wave  processing  than  for  Love  wave  processing.  This  is  due  to  two 
factors.  First,  Love  wave  energy  may  arrive  as  much  as  20°  off  the  great 
circle  azimuth,  and  the  TCA  processor  suppresses  this  motion  as  though  it 
were  noise.  Second,  the  TCA  processor  confuses  the  radial  motion  of  the 
Rayleigh  wave  with  off-azimuth  Love  waves,  and  suppresses  whatever  genuine 
Love  wave  energy  which  is  present.  Thus,  the  TCA  processor  has  shown 
itself  to  be  less  effective  in  separating  Love  waves  than  Rayleigh  waves  from 
noise. 

This  report  describes  an  attempt  to  rectify  this  situation.  The 
TCA  processor  is  modified  to  track  the  incoming  Love  wave  in  azimuth  and 
to  pass  its  entire  duration.  This  is  accomplished  by  using  in  the  design  of 
the  transverse  component  filter  the  additional  information  contained  in  the 
phase  relationships  of  the  Rayleigh  waves. 
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The  evaluation  of  this  modification  of  the  TCA  processor  uses 
both  synthetic  and  real  data.  The  synthetic  data  consists  of  known  signals 
with  high  signal-to-noise  ratios  buried  in  seismic  noise.  These  are  used  to 
study  the  signal-to-noise  ratio  improvement  characteristics  of  the  processor. 
Part  of  the  real  data  used  consists  of  hour-long  noise  samples,  which  are  used 
to  determine  signal  detection  criteria  and  false-alarm  characteristics  of  the 
processor.  The  remainder  of  the  real  data  consists  of  single-site  and  beam 
data  recorded  at  the  Alaskan  Long  Period  Array  (ALPA).  This  data  is  used 
to  study  the  detection  performance  of  the  processor. 

To  allow  direct  comparisons  between  the  original  and  new 
Love  wave  processors,  the  original  Love  wave  processor  was  re-evaluated 
using  the  same  data  base  as  was  used  for  the  evaluation  of  the  new  Love  wave 
processor.  The  Rayleigh  wave  processor  was  also  re-evaluated,  since  the 
necessary  data  processing  for  this  was  performed  automatically  with  the  Love 
wave  data  processing.  Use  of  the  larger  data  base  should  improve  previous 
estimates  of  Rayleigh  wave  detection  capability  improvement  due  to  applica- 
tion of  the  TCA  processor. 

Section  II  of  this  report  describes  the  theory  of  the  TCA  pro- 
cessor, with  emphasis  on  the  new  Love  wave  processor.  The  experimental 
data  base  is  also  described  in  this  section.  Section  IH  describes  the  single- 
site and  beam  processor  gains  to  be  expected  from  the  TCA  processor. 

Section  IV  gives  the  results  of  evaluating  the  detection  capability  improve- 
ment due  to  the  TCA  processor.  Conclusions  regarding  the  method  are  pre- 
sented in  Section  V. 
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SECTION  II 

METHOD  OF  EVALUATION 

A.  DESCRIPTION  OF  THE  THREE-COMPONENT  ADAPTIVE  PROCESSOR 

The  three-component  adaptive  processor  forms  Love  and  Ray- 
leigh wave  filters  in  the  frequency  domain,  which  pass  energy  having  the  phase 
characteristics  of  these  waveforms.  The  filters  are  adaptive,  since  the  data 
are  segmented  and  new  filters  are  designed  for  each  segment.  The  filter 
weights  depend  only  on  the  signal  behavior  during  the  segment. 

The  Rayleigh  wave  filter  algorithm  will  not  be  discussed  in  de- 
tail here,  since  it  is  throughly  described  in  earlier  reports  (Lane,  1973; 
Strauss  and  Tolstoy,  1974).  In  general,  the  method  is  as  follows.  To  pass 
Rayleigh  waves-  the  processor  seeks  signals  which  are  90°  out  of  phase  on 
the  vertical  and  radial  components.  The  Fourier  transform  of  a segment  of 
data  is  taken.  The  phase  angle  0 between  the  radial  and  vertical  components 
is  then  calculated  at  each  frequency.  This  angle  will  be  ff/2  for  pure  Ray- 
leigh motion.  The  filter  weight  at  each  frequency  is  calculated  from: 

F(i)  = sin^  ^ 0 (i) 

(In  this  evaluation,  we  used  N=8.  ) The  Fourier  components  at  each  frequency 
are  then  multiplied  by  the  corresponding  filter  weight.  Finally,  the  data  are 
inverse  transformed.  Processing  is  performed  on  64 -point  segments  with 
fifty  percent  segment  overlap. 

The  computation  method  of  the  original  Love  wave  portion  of 
the  TCA  processor  rests  on  the  fact  that,  when  the  components  are  properly 
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oriented,  the  Love  wave  appears  only  on  the  transverse  component.  The 
Love  wave  filter  weights  are  as  follows: 

G(i)  = cosN|y(i) 

where  y is  the  angle  between  the  transverse  component  and-  the  actual  Love 
wave  particle  motion.  The  exponent  N is  the  same  as  for  the  Rayleigh  wave 
filter.  This  filter  will  pass  in  full  only  that  energy  which  has  the  characteris- 
tics of  on-azimuth  Love  waves. 

The  new  method  of  handling  the  Love  wave  is  to  track  the  in- 
coming energy  in  azimuth  and  to  pass  it  throughout  its  entire  duration.  This 
is  performed  in  the  following  manner  for  each  frequency  component.  First, 
the  time  origin  is  shifted  so  that  the  vertical  component  is  purely  real.  Next, 
the  propagation  direction  (3  (Figure  II- 1)  of  the  transverse  component  energy 
relative  to  the  radial  direction  is  calculated.  Following  this,  the  horizontal 
components  are  rotated  about  the  vertical  axis  by  -p  with  the  result  that  all 
Love  wave  motion  lies  on  the  transverse  component.  After  each  frequency 
component  has  been  rotated  in  this  fashion,  the  time  origin  is  shifted  back. 
Finally,  the  Fourier  transformation  with  fifty  percent  overlapping  of  adjacent 
segments,  as  in  the  original  TCA  processor,  will  yield  the  complete  Love 
wave  time  history.  The  equations  relevant  to  this  process  are  given  in  the 
following  paragraphs. 

After  shifting  the  time  origin  so  that  the  vertical  component  is 
purely  real  (has  zero  phase),  we  have: 


z = 
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P 

(II- 1) 
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(II -2) 
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(II -3) 

R.  = 

i 

L sin 

ib 

sin  p + Pcos  oP  cose 
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Propagation  Direction 
of  the  Rayleigh  Wave 


Propagation  Direction 
of  the  Love  Wave 


Love  Wave  Particle  Motion 


FIGURE  II- 1 


LOVE  AND  RAYLEIGH  ARRIVAL  AZIMUTHS  FOR  AN 
ARBITRARY  FREQUENCY  COMPONENT 


where 


T = L cos  l!i  cos  B 
r 


T.  - L sin  il  cos  p + oPsin  ( oP  sine 


(II -5) 
(II-6) 


Z,  R,  T - refer  to  the  Fourier  components  of  motion  at  the 
same  fixed  frequency  for  the  vertical,  radial,  and 
transverse  components, 

r»  I ■ subscripts  denoting  real  and  imaginary  parts, 

L - Love  wave  amplitude, 

P - vertical  Rayleigh  wave  amplitude, 

o - ellipticity, 

I b - phase  angle  of  the  Love  wave  with  respect  to  the 

vertical  Rayleigh  wave, 

e - propagation  direction  of  the  vertical  Rayleigh  wave 

relative  to  the  radial  diretion, 

p - propagation  direction  of  the  transverse  Love  wave 

relative  to  the  radial  direction. 


Division  of  equation  (II- 3)  by  equation  (II- 5)  yields  the  angle  p: 


r _ Lcos  ll)  sin  p 
Lcos  j|)  cos  p 


= tan  P 


(II-7) 
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Rotation  of  the  horizontal  components  by  the  angle  -P  about  the 
vertical  axis  yields: 
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Under  the  reasonable  assumption  that  € and  p are  about  the  same,  equations 
(11-12)  and  (11-14)  become,  respectively, 
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R.  = aP 

i 

(H-15) 

i 

T.  = L sin  lb  • 

l 

(11-16) 

After  each  frequency  component  has  been  rotated  in  this  way,  the  origin  is 
time  shifted  back  and  an  inverse  Fourier  transformation  with  fifty  percent 
overlapping  of  adjacent  segments  as  in  the  original  TCA  processor  will  yield 
the  complete  Love  wave  time  history.  A listing  and  description  of  usage  of 
the  new  TCA  processor  program  are  given  in  Appendix  A. 

Before  proceeding  to  the  evaluation  of  the  new  Love  wave  TCA 

processor,  two  points  must  be  discussed.  The  first  point  to  be  considered 

is  the  case  where  the  phase  angle  \b  of  the  Love  wjve  relative  to  the  vertical 

Rayleigh  wave  is  ff/2.  Reference  to  equations  (II-3)  and  (II-5)  of  this  section 

shows  that,  in  this  case,  the  angle  p cannot  be  determined,  since  both  T and 

r 

R^  will  be  zero;  therefore  no  rotation  will  take  place.  Although  the  TCA  pro- 
cessed trace  will  not  be  degraded  relative  to  the  unprocessed  trace  in  this 
case,  we  cannot  expect  much  signal -to -noise  ratio  improvement  if  this  is  a 
frequent  occurrence. 
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We  cannot  easily  compute  lb  directly,  but  we  can  check  T 
and  Rr>  Assuming  L is  not  zero,  if  both  Tr  and  Rr  are  zero  simultaneously, 

0 must  be  n/ 2.  Therefore,  the  values  of  Tr  and  R^  as  calculated  by  the  new 
TCA  processor  were  printed  out  for  several  events.  Examination  of  these 
values  did  not  reveal  any  occurrences  of  this  situation.  Thus,  this  matter 
should  pose  no  problem  in  the  operation  of  the  new  TCA  processor. 

The  second  point  to  be  considered  is  that  a high  signal-to-noise 
ratio  on  the  vertical  axis  is  required  to  accurately  pick  the  time  origin  so 
that  the  vertical  Rayleigh  wave  Fourier  component  is  purely  real.  The  Ray- 

« 

leigh  wave  filter  calculates  a measure  of  this  signal-to-noise  ratio  in  the 
phase  angle  between  the  radial  and  vertical  motion.  It  was  originally  pro- 
posed that  when  this  is  near  90°,  little  noise  would  be  present;  therefore,  the 
time  shift  could  be  performed  accurately,  while  in  other  cases  no  time  shifts 
would  occur.  However,  further  consideration  of  this  matter  shows  that  for 
small  signals  (at  or  near  the  noise  level)  the  signal-to-noise  ratio  is  so  low 
that  no  processing  would  occur.  Since  it  is  precisely  these  signals  which  we 
most  want  to  enhance,  no  restriction  of  the  type  described  above  was  used  in 
this  evaluation. 

B.  DATA  BASE 

Suites  of  events  from  the  southern  Kamchatka  and  central  Asia 
seismic  regions  were  selected  for  this  evaluation  of  the  new  TCA  processor. 

Two  regions  were  chosen  rather  than  one  to  allow  a judgment  to  be  made  as 
to  whether  the  effects  of  the  TCA  processor  vary  greatly  from  region  to 
region.  These  regions  were  defined  in  the  final  report  on  ALPA  (Strauss,  1973). 

The  only  restrictions  placed  on  event  selection  were  that  the 
data  were  required  to  be  free  of  spikes  and  long-period  transients  and  unmixed 
with  other  signals  and  that  the  single  site  and  beam  traces  be  4096  seconds  in 
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length.  There  were  a total  of  86  events  available  from  the  southern  Kamchat- 
ka (SKAM)  region  and  76  events  from  the  central  Asia  (CENA)  region.  The 
parameters  of  these  events  are  listed  in  Tables  II- 1 and  II-2,  respectively. 

We  note  that  the  event  numbers  listed  under  the  column  headed  'EVNO'  are 
as  used  by  Lambert  in  his  evaluation  of  the  Very  Long  Period  Experiment  sta- 
tions (Lambert,  et  al.  , 1974).  The  column  headed  'EVENT  NAME'  lists  the 
unique  twelve-character  name  assigned  to  each  seismic  event.  The  columns 
headed  'DATE'  and  'ORIGIN  TIME'  give  the  date  and  time  of  occurrence  of 
each  event.  The  epicentral  coordinates  are  listed  in  the  columns  headed 
'LATITUDE'  and  'LONGITUDE'.  The  column  headed  'DEP'  gives,  when 
known,  the  depth  to  the  hypocenter  of  the  event.  The  column  headed  'MB' 
gives  the  bodywave  magnitude  of  each  event.  The  last  column,  headed  'IS', 
gives  the  information  source  for  the  listed  event  parameters,  where  P stands 
for  the  PDE  (Preliminary  Determination  of  Epicenters)  event  lists,  L stands 
for  the  LAS  A (Large  Aperature  Seismic  Array)  event  lists,  N stands  for  the 
NORSAR  (Norwegian  Seismic  Array)  event  lists,  and  I and  J stand  for  the 
verified  and  unverified  International  Seismic  Month  event  lists,  respectively. 
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SECTION  III 

TC A PROCESSOR  GAINS 
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A.  INTRODUCTION 

In  order  to  determine  TCA  processor  gains,  motion  over  a 
2048  second  segment  containing  a high  signai-to-noise  ratio  ("noise-free")  I 

signal  as  recorded  at  ALPA  was  multiplied  by  a scale  factor  and  added  to  the 
noise  in  a 2048  second  segment  immediately  preceding  the  signal  segment. 

The  TCA  processor  was  then  applied  to  the  resulting  composite  trace.  The 
data  were  also  subjected  to  a bandpass  filter  which  rejected  all  energy  outside 
the  0.  024-0.  059  Hz  (17.  0-41. 5 sec. ) passband.  A wide  range  of  scale  factors 
was  used  for  each  event,  so  that  the  composite  traces  ranged  from  the  case 
where  the  signal  was  completely  buried  in  noise  to  the  case  where  it  was 
clearly  visible.  The  purpose  of  this  procedure  was  to  determine  the  dB  gain 
of  the  TCA  processor  over  the  equivalent  bandpass  filter  as  the  signal  emerges 
from  the  noise.  In  particular,  we  were  able  through  this  procedure  to  make 
estimates  of  the  dB  gain  to  be  expected  when  the  TCA  processor  is  applied  to 
low  signal-to -noise  ratio  events. 

Output  signai-to-noise  ratios  were  calculated  by  dividing  the 
maximum  zero -to -peak  amplitude  found  within  the  arrival  time  window  of  the 
processed  trace  by  the  root-mean-square  (RMS)  value  of  the  noise  in  the  time 
interval  prior  to  the  signal  arrival.  Twenty  times  the  logarithm  to  the  base 
ten  of  this  quantity  is  the  signal-to-noise  ratio: 


SNR 


2°  log  i o 


(Peak  Amplitude  \ 
RMS  Noise  / 


Signal-to-noise  ratios  were  calculated  for  both  TCA  processed  and  bandpass 
filtered  traces.  The  gain  of  the  TCA  processor  over  the  equivalent  bandpass 
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filter  is  the  difference  between  the  signal -to -noise  ratio  of  the  TCA  processor 
and  the  equivalent  bandpass  filter  for  the  same  event. 


Before  proceeding  to  the  determination  of  dB  gain  estimates 
for  the  TCA  processor,  one  further  point  must  be  considered.  In  the  new 
Love  wave  TCA  processor,  some  filtering  of  noise  and  off-azimuth  signals 
can  be  performed  by  requiring  that  the  angle  p be  less  than  some  limit. 

(We  shall  call  this  limit  the  p accept-reject  limit.)  If  a value  of  p greater 
than  this  limit  is  found,  it  will  be  assumed  that  the  motion  is  due  to  noise  or  a 
truly  off-azimuth  signal.  In  this  case,  the  real  components  of  the  horizontal 
motion  and  the  imaginary  components  of  the  transverse  motion  will  be  set  to 
zero.  These  unwanted  components  will  then  not  contribute  to  the  TCA-processed 
motion  in  the  time  domain. 

The  answer  to  the  question  of  what  the  p accept-reject  limit 
should  be  is  not  immediately  obvious.  If  the  limit  is  set  too  low,  off-azimuth 
portions  of  the  Love  wave  (due  to  scattering  and/or  multipathing)  will  be  re- 
jected. If  this  limit  is  set  too  high,  the  desired  filtering  effect  will  be  weak- 
ened. 

An  empirical  solution  to  this  problem  was  found  by  selecting 
several  high  signal -to -noise  ratio  events,  scaling  each  signal  and  adding  it  to 
noise  so  that  the  true  signal-to-noise  ratio  of  the  composite  trace  before  pro- 
cessing was  12  dB,  and  applying  the  new  TCA  processor  repeatedly  while  vary- 
ing the  p accept-reject  limit  from  +_5°  to  +^90°.  The  results  of  this  proced- 
ure are  given  in  Table  III  - 1 . 

We  see  from  this  table  that  the  dB  gain  of  the  Love  wave  TCA 
processor  over  the  bandpass  filter  for  two  of  the  test  events  is  highest  at  5° 
and  for  the  other  two  is  highest  at  10°.  When  plots  of  the  data  produced  by 
this  process  were  reviewed,  it  was  noted  that  Love  wave  energy  after  the  first 
arrival  was  increasingly  suppressed  as  the  p accept-reject  limit  was  de- 
creased. It  was  decided  that  the  best  tradeoff  of  maximum  dB  gain  versus 
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TABLE  III- 1 


dB  GAIN  OF  THE  NEW  LOVE-WAVE  TCA  PROCESSOR 
OVER  THE  EQUIVALENT  BANDPASS  FILTER 
AS  A FUNCTION  OF  THE  p ACCEPT-REJECT  LIMIT 


p Limit 

dB  Gain  of  New  TCA  over  Bandpass  Filter 
for  Test  Events 

Event  172 

Event  505 

Event  693 

Event  1085 

+ 5° 

10.  31 

4.  56 

14.  14 

7.  75 

+ 10° 

10.  46 

9.  06 

7.  67 

6.  92 

+ 15° 

9.  80 

6.  89 

5.  64 

5.  36 

+ 20° 

10.  13 

5.95 

4.  80 

5.  03 

+ 25° 

9.  03 

6.  26 

3.  52 

6.  13 

+ 30° 

7.  13 

5.  30 

1.  81 

6.  06 

+ 35° 

7.  76 

5.  12 

3.  42 

5.  77 

+ 40° 

7.  34 

3.  77 

3.  39 

6.  67 

+ 45° 

6.  92 

3.  69 

3.  14 

7.  41 

+ 50° 

5.  26 

2.98 

3.  12 

6.  93 

+ 55° 

4.  15 

3.03 

3.06 

5.  24 

+ 60° 

3.  93 

2.  92 

3.  09 

5.  01 

+ 65° 

3.  84 

2.  45 

3.03 

1.  85 

+ 70° 

3.  48 

2.  49 

2.  66 

2.  08 

+ 75° 

1.  70 

-0.  17 

2.  36 

2.  39 

+ 80° 

1.  71 

-0.  98 

2.  36 

2.  71 

+ 85° 

0.  84 

-1.  36 

0.04 

1.  92 

+ 90° 

0.  78 

-1. 25 

0.09 

2.  13 

ft 
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suppression  of  later  Love  wave  arrivals  was  at  the  p accept-reject  limit 
o 

of  + 10  . Therefore,  all  further  evaluation  of  the  new  Love  wave  TCA  pro- 
cessor will  use  this  value  for  the  p accept-reject  limit. 

We  note  that,  since  only  ALPA  data  was  used  in  this  evaluation, 
ii  should  not  be  assumed  that  the  above  result  is  applicable  to  data  recorded  at 
other  locations. 


We  shall  now  turn  to  the  determination  of  estimates  of  dB  gain 
due  to  the  TCA  processor.  In  the  following  figures,  which  describe  the 
single-site  and  beam  processor  gains,  the  horizontal  axis  is  the  true  signal- 
to-noise  ratio,  in  dB,  of  the  composite  trace  before  processing.  These  values 
were  obtained  by  multiplying  the  maximum  zero -to -peak  signal  amplitude  by 
the  scale  factor  mentioned  earlier  and  dividing  by  the  RMS  value  of  the  unpro- 
cessed composite  noise.  Twenty  times  the  logarithm  of  this  quantity  is  the 
true  signal-to-noise  ratio. 


True  SNR 


20  log 


Scale  Factor  » Maximum  Amplitude 
RMS  Composite  Noise 


The  vertical  axis  is  the  dB  gain  of  the  TCA  processed  trace  over  the  equiva- 
lent bandpass  filtered  trace  at  the  same  scale  factor. 

We  shall  review  the  evaluation  of  the  original  TCA  processor 
before  considering  the  results  of  the  new  Love  wave  processor.  In  order  to 
have  our  results  directly  comparable  to  those  published  previously  (Lane, 

1973)  we  shall  use  the  same  test  event  (KAM-  199-08AR,  event  number  505) 
as  before  to  evaluate  processor  gains. 

We  note  one  further  point  before  proceeding  to  the  considera- 
tion of  single-site  processor  gains.  In  order  to  facilitate  comparisons  between 
the  old  and  new  versions  of  the  TCA  processor,  it  is  useful  to  have  one  number 
which  represents  the  dB  gain.  We  propose  to  use  for  this  the  dB  gain  at  12  dB 
true  signal-to-noise  ratio,  since  we  are  most  interested  in  the  gain  we  can 


III -4 


expect  from  the  TCA  processor  for  signals  at  or  near  the  noise  level.  (A 
true  SNR  of  12  dB  implies  that  the  maximum  signal  peak  is  at  or  near  the 
amplitude  of  the  maximum  noise  peak.) 

B.  SINGLE-SITE  PROCESSOR  GAINS 

Figure  III  - 1 shows  the  results  of  applying  the  original  TCA 
processor  to  the  single-site  test  event.  All  three  components  follow  the  same 
general  trend.  At  low  true  signal-to-noise  ratios,  where  the  composite  trace 
is  essentially  all  noise,  the  gain  remains  roughly  constant.  As  the  true  signal- 
to-noise  ratio  increases,  the  gain  of  the  TCA  processor  over  the  bandpass  fil- 
ter rises  steadily.  Finally,  at  high  values  of  true  signal-to-noise  ratio,  the 
gain  levels  off.  The  important  point  to  note  on  this  figure  is  that  the  gain  for 
the  Love  wave  processor  does  not  begin  to  increase  until  the  true  signal-to- 
noise  ratio  is  much  higher  than  for  the  Rayleigh  wave  processor  gains.  This 
implies  that  the  original  Love  wave  processor  will  not  be  as  effective  as  the 
Rayleigh  wave  processor  in  detecting  low  signal-to-noise  ratio  events. 

At  12  dB  true  signal-to-noise  ratio,  the  estimates  of  dB  gain  of 
the  original  TCA  processor  are  7.  5 dB  on  the  vertical  component,  6.  5 dB  on 
the  transverse  component,  and  9 dB  on  the  radial  component. 

Figure  III- 2 shows  the  results  of  applying  the  new  Love  wave 
TCA  processor  to  the  single-site  test  event.  The  three  curves  plotted  show 
the  dB  gain  due  to  rotation  about  the  vertical  axis  without  application  of  the  p 
accept-reject  limit  (denoted  by  a dashed  line),  the  dB  gain  due  to  application 
of  the  p accept-reject  limit  without  rotation  about  the  vertical  axis  (denoted 
by  a dotted  line),  and  the  dB  gain  due  to  application  of  the  p accept-reject 
limit  and  rotation  about  the  vertical  axis  (denoted  by  a solid  line). 

This  figure  shows  that  simply  rotating  the  data  results  in  a loss 
in  signal-to-noise  ratio  of  1-2  dB.  The  explanation  for  this  loss  is  that  the 
multipathed  signal  arrives  from  a restricted  range  of  azimuths,  while  noise 
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may  arrive  from  any  azimuth.  Thus,  while  rotation  enhances  the  signal  to  i 

some  degree,  it  more  greatly  enhances  the  noise,  resulting  in  a lower  signal-  t 

to-noise  ratio.  On  the  basis  of  the  above,  one  might  expect  that  the  gain  due 
to  application  of  the  p accept-reject  limit  and  not  rotating  would  be  greater 
than  the  gain  due  to  applying  this  limit  and  rotating.  Figure  III- 2 shows  that 
this  is  not  the  case.  By  applying  the  p accept-reject  limit  and  rotating,  we 
get  a slightly  higher  dB  gain.  The  explanation  for  this  is  that  use  of  the  P 
accept-reject  limit  results  in  the  rejection  of  a large  part  of  the  noise  while 
passing  at  least  part  of  the  signal.  Subsequent  rotation  of  the  passed  frequency 
components  then  results  in  enhancement  of  the  signal  as  was  originally 
intended. 


F rom  the  above,  we  see  that  the  best  form  of  the  new  L^ve 
wave  TCA  processor  appears  to  be  the  combination  of  applying  the  p accept- 
reject  limit  and  rotating  by  -p  about  the  vertical  axis.  Detection  of  the  sig- 
nal using  this  method  is  first  made  at  about  4 dB  true  signal-to-noise  ratio. 

At  12  dB  true  signal-to-noise  ratio,  the  gain  of  this  method  is  9 dB. 

In  summary,  at  low  values  of  true  signal-to-noise  ratio  we  can 
expect  gains  of  about  7-8  dB  for  Rayleigh  waves  and  about  6-7  dB  for  Love 
waves  from  the  original  TCA  processor  operating  on  single-site  data.  In  com- 
parison, we  can  expect  gains  of  about  9 dB  for  Love  waves  from  the  new  TCA 
processor  operating  on  single-site  data. 

C.  BEAM  PROCESSOR  GAINS 

Figure  III- 3 shows  the  results  of  applying  the  original  TCA  pro- 
cessor to  the  beam  traces  of  our  test  event  KAM-199-08AR  (event  number  505). 
The  dB  gains  for  the  vertical  and  radial  components  follow  the  same  general 
trend  as  they  did  in  the  single-site  case,  with  the  difference  that  the  rise  in 
gain  begins  at  higher  values  of  true  signal-to-noise  ratio  for  the  beam  data. 

The  dB  gain  for  the  transverse  component  does  not  show  the  upward  break  of 
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the  single-site  case,  but  gradually  increases  as  the  true  signal -to-noise  ratio 
increases  before  finally  reaching  a constant  value. 

At  12  dB  true  signal-to -noise  ratio,  the  dB  gain  estimates  of 
the  original  TCA  processor  over  the  bandpass  filter  are  0.  5 dB  on  the  vertical 
component,  2.  5 dB  on  the  transverse  component,  and  4.  5 dB  on  the  radial 
component.  These  low  gains,  coupled  with  the  relatively  high  values  of  true 
signal-to-noise  ratio  needed  for  detection  to  occur,  lead  to  the  conclusion 
that  the  original  TCA  processor  should  not  be  as  effective  in  increasing  the 
detection  capability  of  an  array  as  it  is  for  a single  site. 

Figure  III-4  shows  the  results  of  applying  the  new  Love  wave 
TCA  processor  to  beam  data.  Just  as  in  the  single-site  case,  three  versions 
were  applied  to  the  test  event:  one  with  rotation  but  no  p accept-reject  limit, 
one  with  the  p accept-reject  limit  but  no  rotation,  and  one  with  rotation  and 
the  p accept-reject  limit.  The  dB  gains  for  the  version  containing  the  p 
accept-reject  limit  but  no  rotation  are  not  shown  in  Figure  III-4,  since  they 
are  essentially  the  same  as  those  for  the  version  containing  both  rotation  and 
the  p accept-reject  limit.  Although  the  version  containing  only  rotation  out- 
performs the  other  two  versions  at  12  dB  true  signal-to-noise  ratio,  we  note 
that  at  no  point  does  the  gain  of  this  version  exceed  0.  5 dB,  whereas  both  the 
other  two  versions  have  maximum  gains  of  4 dB.  Therefore,  we  judge  the 
versions  containing  the  p accept-reject  limit  to  be  better  than  the  version 
without.  Since  both  versions  containing  the  p accept-reject  limit  yield  essen- 
tially the  same  dB  gains,  we  can  select  either  to  process  the  event  suites.  To 
facilitate  comparison  with  the  single- site  data,  we  will  use  the  version  contain- 
ing both  rotation  and  the  p accept-reject  limit. 

At  12  dB  true  signal-to-noise  ratio,  the  gain  of  the  new  Love 
wave  TCA  processor  is  -1.  5 dB.  However,  we  note  that  the  gain  is  rising 
rapidly  at  12  dB  true  signal-to-noise  ratio,  reaching  approximately  2 dB  at  16 
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dB  true  signal-to-noise  ratio.  Therefore,  in  the  beam  data  case,  we  can 
expect  little  if  any  improvement  in  detection  capability  from  either  the  orig- 
inal or  new  TCA  processors. 

D.  DIFFERENCES  BETWEEN  SINGLE-SITE  AND  BEAM  PROCESSOR 
GAINS 

Table  III  - 2 indicates  that  there  is  a great  difference  between 
the  gains  we  may  expect  when  the  TCA  processor  is  applied  to  single-site 
data  and  when  it  is  applied  to  beam  data.  This  difference  in  performance  im- 
plies some  inherent  difference  between  single-site  and  beam  data  due  to 
beamforming.  This  difference  may  lie  either  in  the  noise,  in  the  signal,  or 
in  both. 

As  recorded  at  an  individual  station,  the  noise  field  consists  of 
non-propagating  noise  and  propagating  noise  with  the  characteristics  of  surface 
waves.  The  beamforming  process  suppresses  the  non-propagating  noise  by  a 
factor  approximately  equal  to  the  square  root  of  the  number  of  stations  in  the 
array.  The  amount  by  which  the  propagating  noise  is  suppressed  is  dependent 
upon  the  arrival  azimuth  of  this  noise;  the  farther  it  is  off  the  beam  azimuth, 
the  more  it  will  be  suppressed.  Thus,  after  beamforming  the  noise  will  have 
a relatively  larger  component  of  on-azimuth  propagating  noise  than  before 
beamforming. 

Now  let  us  consider  how  the  TCA  processor  acts  on  noise.  The 
TCA  processor  searches  for  particular  phase  relationships  among  the  compon- 
ents of  motion.  These  relationships  will  not  be  found  (except  possibly  at  iso- 
lated points)  in  the  non-propagating  noise,  since  in  this  type  of  noise  they  are 
random.  They  also  will  not  be  found  in  the  off-azimuth  propagating  noise, 
since  this  noise  will  not  appear  on  the  proper  components.  The  TCA  processor 
will,  however,  pass  the  on-azimuth  propagating  noise,  since  in  this  case  the 
phase  relationships  will  be  correct. 
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TABLE  III-2 

SUMMARY  OF  GAIN  ESTIMATES 


Gain  in  dB 


Component 

Single  Site 

Beam 

Original  TCA-LR-V 

7.  5 

0.  5 

Original  TCA-LQ-T 

6.5 

2.  5 

Original  TCA-LR-R 

9.0 

4.  5 

New  TCA-LQ-T 

9.0 

-1.  5 

II 


Thus  we  see  that  the  beamforming  process  passes  the  same 
on-azimuth  propagating  noise.  We  cannot,  therefore,  expect  any  appreciable 
gain  when  we  apply  the  TCA  processor  to  beam  data,  since  the  type  of  noise 
it  can  suppress  has  already  been  suppressed  in  beamforming. 


We  shall  now  consider  whether  the  difference  between  single- 
site and  beam  processor  gains  may  be  due  to  some  difference  between  the 
single-site  and  beam  signals.  Before  the  beam  is  formed  from  the  individual 
sites,  differences  in  arrival  time  at  each  of  the  sites  (delay  times)  are  com- 
puted. The  signal  recorded  at  each  site  is  then  time-shifted  by  these  delay 
times  to  time-align  the  signals,  which  are  then  added  and  scaled.  If  the  delay 
times  are  in  error,  due  to  improper  location  of  the  event  epicenter,  the  phase 
relationships  between  the  components  of  motion  will  be  distorted,  resulting 
in  a loss  in  signal  amplitude  when  later  processed  by  the  TCA  processor.  If 
this  is  the  case,  we  should  see  greater  suppression  of  the  beam  signal  than 
the  single-site  amplitude  on  the  TCA  processor  output. 

Table  III  - 3 breaks  down  the  dB  gain  estimates  of  Table  III-2  into 
the  gain  due  to  noise  suppression  and  the  gain  due  to  signal  suppression.  The 
gain  due  to  noise  suppression  is  the  difference  between  the  RMS  noise  in  dB  of 
the  trace  before  TCA  processing  and  the  RMS  noise  in  dB  of  the  trace  after 
processing.  The  gain  due  to  signal  suppression  is  the  ratio,  expressed  in  dB, 
of  the  peak  signal  amplitude  measured  after  TCA  processing  to  the  Deak  signal 
amplitude  measured  before  TCA  processing.  The  net  gain  is  then  simply  the 
sum  of  the  gain  due  to  noise  suppression  and  the  gain  due  to  signal  suppres- 
sion. 


This  table  shows  clearly  that  the  differences  between  single- 
site and  beam  processor  gains  are  due  to  differences  in  the  amount  of  noise 
suppression  and  to  differences  in  signal  suppression. 
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TABLE  III -3 

GAIN  ESTIMATES  AT  12  dB  TRUE  SIGNAL-TO-NOISE 
(FROM  TEST  EVENT  505) 


Component 

Single  s' 

Site  s' 

Beam 

Gain  in  dB 
due  to  Noise 
Suppression 

Gain  in  dB 
due  to  Signal 
Suppression 

Net  Gain 
in  dB 

1 

Original  TCA 

Single  Site 

8.  7 

-1.  2 

7.  5 

LR-V 

Beam 

5.  3 

-4.  8 

0.  5 

Original  TCA 

Single  Site 

10.  2 

-3.7 

6.  5 

LQ-T 

Beam 

6.4 

-3.9 

2.  5 

New  TCA 

Single  Site 

12.  2 

-3.  2 

9.0 

LQ-T 

Beam 

4.  7 

-6.  2 

-1.  5 

Original  TCA 

Single  Site 

11.2 

-2.  2 

9.0 

LR-R 

Beam 

7.  9 

-3.4 

4.  5 

SECTION  IV 

DETECTION  CAPABILITY  ESTIMATES 


A.  INTRODUCTION 

Before  processing  the  event  suites  with  the  new  Love  wave 
TCA  processor,  it  was  necessary  to  check  our  definitions  of  detection  crite- 
ria. Therefore,  we  processed  twenty  noise  samples  with  this  method.  The 
output  showed  that  TCA  processed  noise  tends  to  be  predominantly  18-20  sec- 
ond period  energy  with  packets  of  lower  frequency  energy  appearing  as  the  P 
accept-reject  limit  increases.  Signal-like  packets  were  occasionally  found; 
however,  these  did  not  display  dispersion  and  changed  only  slightly  as  the  P 
accept-reject  limit  was  increased.  In  contrast  to  this,  it  was  noticed  that 
signals  appear  to  'grow'  as  the  P accept-reject  limit  increases,  indicating 
the  inclusion  of  more  off-azimuth  (multipath)  energy. 

From  the  above,  it  is  obvious  that  the  principal  detection  cri- 
terion remains  the  presence  of  dispersion  in  the  signal  gate.  Other  detection 
criteria  are  that  the  suspected  signal  begins  close  to  its  predicted  arrival 
time  and  that  the  largest  peak  of  the  suspected  signal  be  3 dB  or  more  above 
any  other  peak  in  a gate  starting  100  seconds  before  the  estimated  signal 
arrival  time  and  ending  100  seconds  after  the  estimated  signal  end  time. 

(This  last  criterion  should  not  be  considered  to  be  absolute.  It  is  occasion- 
ally possible  to  detect  a signal  from  its  dispersion  characteristics  even 
though  its  largest  peak  amplitude  is  equal  to  that  of  the  largest  noise  peaks.) 

All  events  listed  in  Tables  11-1  and  II - 2 were  processed  by  the 
new  Love  wave  TCA  processor  and  evaluated  with  the  above  detection  criteria 
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to  determine  the  improvement  in  detection  capability  due  to  this  method.  In 
addition,  since  only  a small  suite  of  events  was  used  in  the  original  TCA  eval- 
uation (Lane,  1973)  and  since  we  wish  to  make  direct  comparisons  between  the 
original  and  new  Love  wave  TCA  processors,  both  event  suites  were  also  pro- 
cessed by  the  original  method.  The  version  of  the  new  Love  wave  TCA  pro- 
cessor used  here  is  the  one  with  rotation  and  the  p accept-reject  limit,  since 
in  Section  III  we  found  that  this  version  yielded  higher  dB  gains  than  the  version 
which  contained  only  rotation  or  the  version  which  contained  only  the  p accept- 
reject  limit. 


B.  SINGLE-SITE  DATA  RESULTS 

We  will  consider  detection  capability  improvement  for  single- 
site data  in  terms  of  the  fifty  percent  detection  threshold  and  in  terms  of  the 
change  in  the  number  of  events  detected.  The  fifty  percent  detection  thresh- 
old is  defined  as  the  value  at  which  the  probability  of  detecting  an  event  is 
fifty  percent.  We  will  not  discuss  detection  capability  improvement  for  the 
radial  component,  since  results  for  this  component  are  essentially  identical 
to  those  for  the  vertical  component.  Figures  B-l  to  B-5  of  Appendix  B pre- 
sent the  detection  statistics  for  the  86  events  of  the  southern  Kamchatka  suite. 
Figures  B-6  to  B-10  present  the  detection  statistics  for  the  76  events  of  the 
central  Asia  suite.  These  detection  statistics  are  summarized  in  Table  IV-1. 

From  Table  IV-1,  we  see  that  the  Rayleigh  wave  part  of  the 

original  TCA  processor  yielded  a fifty  percent  detection  threshold  improvement 

of  0.41  m,  units  for  the  southern  Kamchatka  event  suite  and  0.24  m,  units  for 
b b 

the  central  Asia  event  suite.  These  correspond  to  processor  gains  of  approxi- 
mately 8 dB  and  5 dB,  respectively.  These  values  compare  fairly  well  with  the 
7.  5 dB  gain  estimate  made  in  Section  III.  The  Love  wave  part  of  the  original 
TCA  processor  yielded  a fifty  percent  detection  threshold  improvement  of  0.  17 

m units  for  the  southern  Kamchatka  event  suite  and  0.25  m,  units  for  the 
b b 
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TABLE  IV- 1 

SUMMARY  OF  50  PERCENT  DETECTION  THRESHOLD  ESTIMATES 

IN  UNITS 
b 


Type  of  Estimate 

. 1 

Region 

Southern 

Kamchatka 

Central 

Asia 

LR-V  Bandpass 

4.  61 

4.  59 

LQ-T  Bandpass 

4.  49 

4.73 

LR-V  Original  TCA 

4.  20 

4.  35 

LQ-T  Original  TCA 

4.  32 

4.  48 

LQ-T  New  TCA 

4.  11 

4.  38 
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central  Asia  event  suite.  These  correspond  to  processor  gains  of  approxi- 
mately 3.  5 dB  and  5 dB,  respectively.  These  values  compare  fairly  well  with 
the  6.  5 dB  gain  estimate  made  in  Section  III. 

Table  IV-1  shows  that  the  new  Love  wave  TCA  processor  (rota- 
tion and  (3  accept-r eject  limit)  yielded  a fifty  percent  detection  threshold  im- 
provement of  about  0.  35  mfa  units  for  both  regions,  corresponding  to  a pro- 
cessor gain  of  7 dB.  This  compares  quite  well  with  the  9 dB  gain  estimate 
made  in  Section  III. 

Another  way  of  viewing  the  detection  capability  improvement  of 
the  TCA  processor  is  to  consider  the  number  of  events  detected  on  the  pro- 
cessor output  which  were  not  detected  on  the  bandpass  filter  output.  These 
results  are  summarized  in  Table  IV-2.  From  this  data  we  see  that^ne  Ray- 
leigh wave  part  of  the  original  TCA  processor  and  the  new  Love  part  of  the 
TCA  processor  perform  somewhat  better  on  the  southern  Kamchatka  data  than 
on  the  central  Asia  data.  The  original  Love  wave  part  of  the  TCA  processor 
performs  equally  well  on  data  from  the  two  regions.  In  agreement  with  previ- 
ous results,  we  see  that  the  new  Love  wave  part  of  the  TCA  processor  is 
superior  in  performance  to  the  original  Love  wave  part. 

The  data  listed  in  Table  IV-2  indicate  that  both  the  original 
Rayleigh  wave  processor  and  the  new  Love  wave  processor  can  be  expected  to 
detect  approximately  20  to  35  percent  of  those  events  not  detected  on  the  band- 
pass filter  output.  This,  of  course,  assumes  that  the  original  data  base  has 
m^  and  detection/non-detection  distributions  similar  to  the  data  bases  used  in 
this  evaluation. 

C.  BEAM  DATA  RESULTS 

The  results  reported  by  Lane  (1973)  and  the  results  of  Section 
III  of  this  report  all  indicate  that  the  TCA  processor  in  both  its  original  and 
new  forms  is  of  little  value  when  applied  to  beam  data.  For  the  sake  of 
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thoroughness,  however,  the  beams  of  the  central  Asia  suite  of  events  were 
processed  in  the  same  manner  as  the  single-site  data.  The  results  of  this 
processing  are  shown  in  Figures  B-ll  to  B-15  of  Appendix  B and  summarized 
in  Table  IV-3.  These  results  simply  reconfirm  previous  results.  The  only 
significant  improvement  in  detection  capability  was  made  on  the  LR-V  part  of 
the  TCA  processor.  This  0.  15  m^  unit  lowering  of  the  fifty  percent  detection 
threshold,  in  conjunction  with  the  insignificant  changes  in  detection  threshold 
due  to  either  form  of  the  Love  wave  TCA  processor,  indicate  that  it  would  not 
be  productive  to  routinely  apply  the  TCA  processor  to  beam  data. 


SECTION  V 


CONCLUSIONS 


The  following  results  have  been  found  during  the  course  of  this 
evaluation  of  the  Three -Component  Adaptive  (TCA)  processor: 

• The  most  effective  form  of  the  new  Love  wave  TCA  processor 
contains  rotation  of  the  Love  wave  frequency  components  about 
the  vertical  axis  to  the  transverse  component  with  an  accept- 
reject  limit  placed  on  the  Love  wave  arrival  azimuth. 

• At  12  dB  true  signal-to-noise  ratio,  we  can  expect  7-9  dB  gain 
for  Rayleigh  waves  and  6-7  dB  gain  for  Love  waves  from  the 
original  TCA  processor  on  single-site  data. 

• At  12  dB  true  signal -to -noise  ratio,  we  can  expect  about 
9 dB  gain  for  Love  waves  from  the  new  TCA  processor  on 
single-site  data. 

• Only  low  or  negative  gains  can  be  expected  from  either  version 
of  the  TCA  processor  when  applied  to  beam  data. 

• Use  of  the  original  TCA  processor  on  single-site  data  lowered 
the  fifty  percent  Rayleigh  wave  detection  threshold  by  about 

0.  25-0.  4 units  and  the  fifty  percent  Love  wave  detection 

threshold  by  about  0.  15  m,  units. 

b 

• Use  of  the  new  Love  wave  TCA  processor  on  single-site  data 

lowered  the  fifty  percent  Love  wave  detection  threshold  by  about 

0.  35  m,  units, 
b 
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• Use  of  the  original  TCA  processor  on  beam  data  yielded  a 

small  (0.15  m unit)  improvement  in  the  Rayleigh  wave  fifty 
b 

percent  detection  threshold  and  essentially  no  change  in  the 
Love  wave  fifty  percent  detection  threshold. 

• Use  of  the  new  Love  wave  TCA  processor  on  beam  data  yielded 
essentially  no  change  in  the  Love  wave  fifty  percent  detection 
threshold. 

The  above  points  make  it  clear  that  the  original  Rayleigh  wave 
TCA  processor  and  the  new  Love  wave  TCA  processor  are  effective  when 
applied  to  single-site  data.  No  version  is  effective  when  applied  to  beam  data. 
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APPENDIX  A 

THE  UPGRADE  TCA  PROGRAM 


The  program  listing  given  on  the  following  pages  is  the  most 
general  version  available.  It  will  work  on  either  single-site  or  beam  data 
and  can  be  used  to  determine  processor  gains  or  to  enhance  detection  studies. 
To  use  this  program,  note  the  following  data  card  inputs: 

DATA  CARD  INPUT 


Card  Column  Format 


Parameter 


1-12 

3A4 

Event  name  to  be  processed 

21-22 

12 

Vertical  Trace  Number 

23-24 

12 

Transverse  Trace  Number 

25-26 

12 

Radial  Trace  Number 

5 

11 

ISUM.  If  ISUM=0,  process  entire 
trace.  If  ISUM^O,  scale  second 
half  of  trace  and  add  to  first  half. 

6-15 

F10.  3 

SCFAC  - factor  with  which  to  scale 
second  half  of  trace  (used  if  ISUM 

t 0). 

16-25 

F10.  4 

FLO  - low  frequency  of  bandpass 
filter 

26-35 

F10.4 

FHI  - high  frequency  of  bandpass 
filter 

36-45 

F10.  2 

ANGBET  - accept/ reject  limit  for 
Love  wave  arrival  azimuth 

50 

11 

IFLAG.  IfIFLAG=0,  process 
single-site  data;  if  IFLAGi/0,  pro- 
cess beam  data. 

In  general,  the  trace  numbers  used  on  card  one  arc  as  follows: 


Trace  1 - vertical  single-site  data 
Trace  2 - transverse  single-site  data 
Trace  3 - radial  single  site  data 
Trace  4 - vertical  beam  data  (NORSAR); 

transverse  beam  data  (ALPA)  1 

Trace  5 - transverse  beam  data  (NORSAR); 

vertical  beam  data  (ALPA)  1 

Trace  6 - radial  beam  data 

Input  Tape  - TDFILTER  beam  tape;  mounted  on  SYSG18  1 

Output  Plot  Tape  - Scratch  tape  mounted  on  SYS006  ® 

Restrictions: 

All  three  components  must  be  available.  I 

Data  length  must  be  a power  of  two.  _ 

I 
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MAIN  PROGRAM 
(PAGE  1 OF  2) 

DIMENSION  £HEAD(400) ,5EEM(210e,3) , IEHEAD(400) 

OlhtfcS  tvJN  EVENTC3) ,THEAD(40) 

01  •■'E  MS  I ON  DTYPt  (2) 

DIMENSION  c O I ( 2 0 ) 

I NTEGF  R *2  TAPENUC3) 

INTEGER  BEAMN0C3) 

EQUIVALENCE  (TAPEnO, THEAD(2)) 

EQUIVALENCE  (EHEAD, IEHEAD) 

DATA  OTYPE/'TOBE', 'AM  '/ 

U AT  A EOl/'EOItf', 19*'  '/ 

ITRACRC  J)=3*(  J-D+2 
10  = 1 8 

CAL'.  TIOISK('OPEN') 

CALL  TI0PE\(IU,THEAD,RP99P) 

**  I TEC  3,  SO  TnEAD 

50  Fjr<MAT(/////2PX,  'INPUT  TAPE  LABELS  ARE  ' // , 2 ( 1 0X , 1 H ' , 20 A4, 1 H' , / ) ) 

RE  AC/  C HTPOL  CAROS 

100  REAu(t,  U«,£N0s9998)  EVENT,  (BEAMNO(I),  1=1,3) 

110  F0k:UT(3AU,9X,3I2) 
hR1TEC3, 120)  EVENT 

120  FORMAT  f 1H1,  10X,  'EVENT  NAME  = ' , 1 H'  , 3 A4,  1 H'  ) 

CALL  T I INTPC' INPUT'. 'SEARCH', IU,IEHEAD, EVENT,  IC) 

IF (IC,E*.0)  GO  TO  130 
WR I TE ( 3 , 126)  IC 

126  FORMAT (10X, 'TIIhTP  ERROR  CODE  = ',Z8) 
wR ITE ( 3, 125)  EVENT 

125  FORMAT ( 1 h 1 , 1HX, 'EVE’ iT  ',1H',3A4,'  NOT  FOUND--ABORT' ) 

GO  TO  9998 

130  CALL  DTSrI0('«RIT',1, IErtEAD, 340, 1100) 

wRITtCi, 121)  BEAMNO(1),BEAMNO(2),0EAMNO(3) 

121  FORMAT ( 1 0 X , 'BEAMS  TO  PROCESS  (V,T,R)  = ',312) 

C*lL  L S H E A J ( 4 , I E H £ A 0 ) 

NBEEMSr IEHE AD(29) 

TIMES=IErtEA0(46) 

NeyTsIEHEdDCB) 

L£N=(NbYT/4)-25 
TlMEL  = IEMi:Ai)(6B)-IEHtAD(46) 

TIMt.P  = IfcrtEAD(70)-IEHEAD(46) 

TIMEE=T1MER+IEHEAD(73) 





MAIN  PROGRAM 
(PAGE  2 OF  2) 

C 

C READ  INPUT  DATA  FROM  TAPE  AND  PUT  ON  DISK 

C 

P'J  220  I*  1 » NBEtMS 

CALL  T I*E  AD  ( I'J,  BEEM,  NBYT,  IC,  &16F,  8 1«0,&160) 

CALu  0I3KI0C  'wR  IT',  1 TRACK  (1),  BEE  M,NBYT/4, 8,200) 

GO  TO  220 

160  «RITt(3,  1 7 y ) I,  IC 

170  FORMAT (/2X, * T APE  ERROR  ON  BEAM  NUMBER  s' , IS, 5X , ' IC= ' , 28) 

GO  TO  220 

ISO  RRlTfe(j,lR,i)  I 

190  FOR-.AT  (//'iX, 'aEAM  NUMBER  s', 13,'  NOT  FOUND') 

GO  TO  220 

200  *RITE(3, 2H’)  1 

210  FOR’  AT (//5X, 'DISK  ERROR  On  wRITE  INDEX  s', 13,//) 

222  CONTINUE 
C 

Du  2 7'*  n s 1 , 3 
Do  2 3 ‘‘  1*1,2100 
23i  BEE  m ( T , >\)  =o , y 

NijMbEmsBEAMNO(K) 

IF(NU*‘RfcM.UE.NBEEMS)  GO  TO  290 
R 1 1 1 ( 3 , 1 R 0 ) NU'IUEM 
GO  TO  270 

29/  CALL  DISK  IOC 'READ',  I TR A CK ( NUMBEN) , BEEM  ( 1 , K ) , NB  Y T/ 9 , & 2E  v; 

GO  TO  270 

| 252  WRITE  (3,260)  K,nU'1B£M 

26/  FOR-' ATC//5X, 'DISK  ERROR  ON  READ.  LOOPS', 13,'  NUMBEMs  ' , T V/ ) 
270  CONTI  <JC 

CALL  LMO A Dl( BEEM, lEN, EVENT, TIMES, TIMED, TIMER, TIME  " 

GO  TO  1-5  2 

999«  CALL  T TOISK('CLQS') 

CALL  PLOT (4. 0,0. 0,999) 

R999  CALL  T ITAME('UNLOAD',lU,BEEM,NBYTE, IC) 

STOP 

END 


! 
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SUBROUTINE  LHQADL 
(PAGE  1 OF  5) 

SUBROUTINE  LHCA0L(DATA,-!TCIM,NAME,T1MES,7IMEL,TIMER,TIMEE) 
DIMENSION  F (3),TEmP(3),!>ATA(2100,3),PROC(2100,3),IBIT(512) 
DIMENSION  T RG ( 5 1 2 ) » 1 1 B I T ( 5 1 2 ) , TP IG (51 2) , PBUF(  64,  3,  2) 

DIMENSION  FIN(2106*,3) , TSQO ( 3) , CT IME ( 3) , NAME ( 1 ) , ASK (20) 

OAT  A ASK/265*'****'/ 

IPTSEGsbu 
IFO^E  f<  — 3 

100  CALL  PICT  (f .0,-12. 0,-3) 

CALL  PLOT  (Q.0, ltf.0,-3) 

CALI  PLOT  0>. ,-0.5, -3) 

ISOM  EG  0 PROCESS  ENTIRE  TRACE 

I SON  NE  65  SCALE  SFCOND  HALF  OF  TRACE  AND  ADD  TO  FIRST  HALF 

SCFAC  IS  the  SCALE  FACTOR  FOR  THE  SIGNAL  HALF  OF  THE  TRACE 

IF  I S U M.  = 3,  SCFAC  MUST  = 1.0 

FLO  = L OH  BANDPASS  FREQUENCY 

FHI  = HIGH  BANDPASS  FREQUENCY 

ANoBFT  = ACCEPT/PEJECT  LIMIT  ON  BETA 

IFlAG  EG  0 SINGLE  SITE  DATA 

IFLAG  NE  0 PEAR  DATA 

TAPR  r AMOUNT  OF  BP  FILTER  TAPER--MAY  NOT  BE  ZERO 

REACT  1, Ilf, ENDSP9R0)  ISUM,SCFAC,FLO,FHI,ANGBET,IFLAG,TAPR 
110  FORMAT  CIS,  FI  C.  3#  Fl  0. 4,  F 16).  4,  Fl  0.2,  IS,  FI  0,4) 

CALL  SY*BOL(0.0*-2,0#fc.2,'FVENT  NAME', 0.0, 10) 

CALL  SYMBOL  CO. 65, -3.6 :,a,2,  'START  TIME'.  0. 0,  10) 

CALL  SYMBOLC0.0,-«.0,F.i5,'LON  BP  FREQ', 0.0, 11) 

CALL  SYMBOLS. 0, -4.5,6'. 15,  'HIGH  BP  FREQ' , 0, 0,  1 2) 

CALL  MHBER(2.5,-4,5,C?,  ’ 5,  F HI,  0. 0, 4) 

CALL  NUMBER (2. 5, -4, 0,0.  5, FLO, 0,0, 4) 

CALL  NUMBER (2. 5, -3. 0,0.2, TIMES, 0, 0,-1) 

CALL  SY'  6CL(2.5,-2.0,6'.2,NAME,0.0,  12) 

IFITFLAG.NE.P)  GO  to  115 

CALL  SYMBOL  CO. P, -5.0, 0, 15, 'SINGLE  SITE  DATA', 0, 0, lfe) 

GO  TC  117 

115  CALL  SYHBOLCf*. 0,-5, 0,0.  15,  'BEAM  DATA', 0,0, 9) 

117  IF  (ISUM.EG.k5)  GO  TO  120 

CALL  SYMbOLCtf. 0,-5. 5,0. 15, 'SIGNAL  S, F. ', 0.0,  1 1 ) 

CALL  M:MBER(2.5,-5.5,P. I 5, SCFAC, 0. 0, 4) 

120  CALL  PLOT (12. 0,1. 4, -3) 

SCALE  SIGNAL  AND  ADD  TO  PRECEEDING  NOISE,  IF  DESIREO 

IFUSUf  .EO.0)  GO  TO  215 
TINELs , IMEL-NTOIH 
TIHEkstIMER-NTDIm 
TINEErTlMEE-NTOlK 

CALL  SC.  ALA  0 (CAT  A , PKOC , K TD I M,  SCFAC,  TIMED 
GO  TO  214 
215  DO  210  J = 1 , 3 

DO  210  I=1,NTD1M 
210  PROCCI, J)=DATA(I,  J) 
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SUBROUTINE  LHQADL 
(PAGE  2 OF  5) 

C 

C COMPUTE  am;  print  out  composite  s/n  ratios 

c 

2)4  CALL  Rr'3(FkOC*CTl;E,TS0D,TIMEL,TlMEE,F,TEMP,ISUM) 
*K17fc(3*216HAS*d>»  1=1*2/) 

216  F0M'ATC1S*,?FAU) 

«MTF(?.21I)  FLO*  FHI * TAPR 

211  FUR’ AT ( iva, '*****  FLO  = ',F7.4,'  **  FrtI  = ',F7.4,'  **  TAPER  = 

1 ' , F 7 . 4 , ' ***♦*') 

IF  ( I b 0 Kl , t < 1 • / ) GO  TU  219 
k*ITE(3, 2) 7) 

217  FOP ’ AT ( 5X , ' ') 

wRITE(3*?l&)  SCFAC 

218  FQRMATC35X, 'SIGNAL  SCALING  FACTOR  = ',F7.b) 

* R I T E ( 3 * 2 1 7 ) 

219  *HJTE(3*22iJ) 

220  FCFmaU2/a, 'UNPROCESSED  DaTA') 

Kr  l T E ( i * 2 3 »>  j 

23k*  FORMAT  (?•/*,  'S/N',  12X,  'AT  T I ME  ' , 6X  , 'RMS  NOISE'*  1 H X , 'PEAK') 

WRITE  (3#?**0)  (TSL’DC  J)  , CTIMfc(J)  ,F(  JJ  , TEMP C J J , J= 1 , 3 ) 

242  FORMAT  (2C  X,  Ffc.2*9Xf  F6.W,  U'X*F5.2*  1 HX* E 1 0. 3/ * 20X , F6. 2, 9X, F6 . F , 10X,F 
15.2* 1FX*F_ if .3/*20X*F6.2,9X,F6.0, IBX,F5.2* 1<3X,E10.3) 
kRJTF (3*21  7) 

PERFORM  THREE-COMPONENT  ADAPTIVE  PROCESSING 

DO  2 Of.  J = 1 * 3 
DO  2b'  1=1, 1PTSEG 
25'J  PFUFCI,  J,  l)=PROC(I,  J) 

262  CALL  T J OFF  7 (PdUF (1*J* 1 ) * TRG* IB1T*  IPTSEG) 

CALL  LQVFIUPRUF,  1PTSEG,  ANGBE7) 

CALL  PATFILCPBUF* IPTSEG* IPOwEH) 

DO  33/  .1  = 1*3 

332  CALL  TIIFF  7 (PBUF(1 *J* 1),TRIG» IIBIT* IPTSEG) 

Iu2=IPTSF G/2 
D 1 1 3 a c .'  .1  = 1*3 
DO  340  1=1,102 
342  F I M I * J ) sPntlp  (I,  J,  1) 

1 UT C C.  = ( MO  i M/ IP  7SEG  1*2-2 
NCCnf TsIU2 
DO  abf  kk.KsI  , INTOO 
DO  350  J = 1 * 3 
DO  3b1/  1 = 1*  IPTSEG 
PHLFII* J,2)=P0UF(  I, J, 1) 

342  PhuKl,  J,  1 ) =PRQC ( I + nCC’UNT  , J) 

DC  36^  j = 1 * 3 

360  CALL  TIlFFTCPBUFd,  J,  1)*TRG.  IBIT, IPTSEG) 

CALL  L&vFILCPSUF* IPTSEG* AN&BET) 

CALL  RAVFILCPftUF* IPTSEG, IFOwER) 

DO  43/  J=l*i 

4 3C5  CALL  TlIFFT(PHuF(l,J,l)*7RI&*IIBIT, IPTSEG) 
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EI02*I02 
DO  440  1*1 , 102 
ECO'JNTsI 
RAT*EC0UNT/EI02 
DO  440  J*  1 , 3 

440  FINCItNCOUNT, J)*PBUF(I, J, 1 ) *R AT+PBUF ( I+I02, J, 2) * ( 1 , 0-RAT) 

450  NC0UNT*NC0UNTtI02 
C 

C BANDPASS  FILTER  THE  ADAPTIVE-PR0CE88ED  TRACES 
C 

DO  460  J*  1 , 3 
DO  460  1*1,102 

460  FINCI+NCOUNT, J)SPBUF(I*X02, J, 1) 

DO  470  J«i , 3 

470  CALL  TIDFFT(FIN(1 , J) , TR6, IBIT,NTDIM) 

CALL  FILTBP(FIN,NTDIM,FLO,PHI,TAPR) 

DO  500  J*1 , 3 

500  CALL  TIIFFT(FIN(1,J) , TRIG, I1BIT,NTD1M) 

COMPUTE  AND  PRINT  OUT  ADAPTIVE  8/N  RATIOS 

CALL  RMS (FIN, CTIME,TSQD,TIMEL,TIMEE,P, TEMP, ISUM) 

WRITE(3,510) 

510  FORMAT (20X,#TCA  PROCESSED  DATA*) 

WRITE (3, 230) 

WRITE(3,240)(TSQD(J),CTZME(J) , F( J) , TEMP ( J) , J«1 , 3) 

WRITE (3, 515)  ANGBET 

515  FORMAT (15X, ******  LQ-T  ACCEPT/REJECT  LIMIT  ON  BETA  ■ ',F6,2,*  DEGR 
1EES  ***•*') 

PLOT  TCA  PROCESSED  TRACES 

CALL  SYMBOL (-3,0, -1,0, 0,2, 'VERTICAL  ADAPT*, 0,0, 18) 

CALL  SYMBOL(-3.0,-1.6,0.15,'S/N',0,0,3) 

CALL  NUMBER (-2, 0,-1, 6,0, 1 5, TSQD ( 1 ) , 0, 0, 2) 

CALL  SYMB0L(-3, 0, -1,9,0, 15, *AT  TIME*, 0,0, 7) 

CALL  NUMBER (-1 ,0,-1, 9,0, 15, CTIME(1),0,0,-1) 

CALL  PLOT (0,0, 0,0, 3) 

CALL  PLTRTN(FIN(1,1),NTDIM,2,0,TIMEL,TIMER,TIMEE,SF) 

CALL  SYMBOL(-3.0,0.3,0,15,*3.F,*,0,0,4) 

CALL  NUMBER (-2, 4, 0, 3, 0, 15, SF, 0,0, 6) 

CALL  PLOT (0,0, 0,0, 3) 

CALL  SYMBOLC-3, 0,-1. 0,0. 2, 'TRANSVERSE  ADAPT', 0,0, 16) 

CALL  SYMBOL(-3,0,-1.6,0.15,'8/N',0,0,3) 

CALL  NUMBER(-2.0,-l,6,0.15,T3QD(2),0.0,2) 

CALL  SYMBOLC-3.0,-1.9,0.15,*AT  TIME', 0,0, 7) 

CALL  NUMBER (-1 , 0,-1. 9,0. 15, CTIME (2), 0. 0,-1) 

CALL  PLOT (0,0, 0.0, 3) 

CALL  PLTRTN(FIN(1,2),NTDIM,2,0,TIMEL,TIMER,TIMEE,SF) 

CALL  SYMBOL(-3,0,0,3,0.15,*$,F,*,0.0,4) 
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CALL  NUMBER (-2, 4, 0, 3,0, 15, 8F, 0,0, 6) 

CALL  PLOT(0.0,0.0,3) 

CALL  3YMBOL(1.0,0.2,0,1, 'BETA  ACCEPT-RE JECT  LIMIT  EQUALS', 0.1, 11) 
CALL  NUMBER(4.2,0,2,0.1,ANGBET,0.0,2) 

CALL  PLOT(0,0,0.0,3) 

CALL  3YMBOL(-3.0,-l«0»0.2, 'RADIAL  ADAPT', 0,0, 12) 

CALL  3YMBOL(-3.0,-t.6,0.15,'8/N',0,0,3) 

CALL  NUMBER(-2.0,-1.6,0,15,T8QD(I),0,0,2) 

CALL  SYMBOL(-3.0,-1.9,0.15,'AT  TIME*, 0,0, 7) 

CALL  NUMBER (-1.0,-1. 9,0. 15, CTIMEC3), 0.0, -1) 

CALL  PLOT(0,0,0,0,3) 

CALL  PLTRTNCFINCl, 3), NTDIM, 2,0, TIMCL, TIMER, TIMES, 8P) 

CALL  SYMBOL(-S,0,0,3,0.t5,'8,F,',0,0,«) 

CALL  NUMBER (*2,4,0, 3,0, 15, SF, 0,0, 8) 

CALL  PLOT (0,0, 0,0, 3) 

BANDPASS  FILTER  THE  COMPOSITE  TRACES 

DO  520  J*l,3 
DO  520  I»1 , NTDIM 
520  FIN(I,J)«PROC(I,J) 

DO  530  J«l,3 

530  CALL  TIDFFT(PIN(I, J) ,TRG, IB IT, NTDIM) 

CALL  FILTBP (FIN, NTDIM, FLO, FHI,TAPR) 

DO  560  J«l,3 

560  CALL  TIIFFT(FIN(1,J), TRIG, IIBIT, NTDIM) 

COMPUTE  AND  PRINT  OUT  BANDPASS  S/N  RATIOS 

CALL  RMS (FIN, CTIME,TSQD,TIMEL,TIMEE,F, TEMP, I8UM) 

MRITE(3,570) 

570  FORMAT(20X, 'BANDPASS  FILTERED  DATA') 

WRITE (3, 230) 

WRITE (3, 240) (TSQD(J),CTIME(J),F(J),TEMP(J),J«1,3) 

CALL  SYMBOL(-3,0, -1.0, 0,2, 'VERTICAL  BP  FILT', 0,0, 16) 

CALL  3YMBOL(-3.0,-l,6,0.15,'8/N',0,0,3) 

CALL  NUMBER (-2,0, •1,6,0, 15, TSQD(l) , 0,0, 2) 

CALL  PLOT(0,0, 0,0, 3) 

CALL  PLTRTN(FIN( 1,1), NTDIM, 2, 0,TIMEL, TIMER, TIMEE,SF) 

CALL  SYMBOL(-3,0,0,3,0.15,'8,F,',0,0,4) 

CALL  NUMBER (-2. 5, 0, 3,0, 15, 8F, 0,0, 6) 

CALL  PLOT (0, 0, 0« 0, 3) 

CALL  SYMBOL(-3.0, -1.0, 0.2, 'TRANSVERSE  BP  FILT', 0,0, IS) 

CALL  SYMBOL (-3,0, -1.6, 0.15, 'S/N', 0,0, 3) 

CALL  NUMBER (-2. 0,-1 .6, 0. 1S,T8QD(2) ,0,0,2) 

CALL  PLOT (0.0, 0.0, 3) 

CALL  PLTRTN(FIN(1,2),NTDIM,2,0,TIMEL,TIMER,TIMEE,8F) 

CALL  SYMBOL(-3.0,0,3,0.15,'S,F.',0,0,4) 

CALL  NUMBER (-2. 5, 0. 3,0, 15, SF, 0,0, 6) 

CALL  PLOT(0,0,0.0,3) 
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CALL  SYMBOL(-3.0r-1.0'0*2»#RADIAL  BP  FXLT',0,0, 14) 
CALL  SYM0OLC-3.0,-1.6,0.l5»#3/N',0.0,3> 

CALL  NUMBER(-2.0,-1.6,0.l5»T3QD(3)#0,0#2> 

CALL  PLOTC0.0,0.0,3) 

CALL  PLTRTN(FIN(1,3),NTDIM,2,0,TXMEL,TIMER,TXMEE,8F) 
CALL  SYMBOL(-3,0,0.3,0,l5,'S,F,',0,0,<») 

CALL  NUMBERC-2,5,0.3,0,15,3F,0,0,6) 

CALL  PLOTC0V0,0.0,3) 

*****  TERMINATE  ***** 

ENTDIMbNTDXM 
FINUP*(ENTD I M/50,0) *4,0 
CALL  PLOT (FXNUP»  0,0 »»3) 

9998  CALL  PLOT  (0,0,0V0*990) 

RETURN 

END 


! 

f I 
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SUBROUTINE  3CALAD(DATA,PR0C,NTDIM,3CFAC,TILQ) 
DIMENSION  DATA(2100,3),PROC(2100,3) 

DIMENSION  TEMPC2100,3) 

ISTART.TILQ/2, 

NTOIMbNTOIM/2 
DO  10  J«l,3 
00  10  I«1,NTDIM 
TEMP(I, J)«0,0 
10  CONTINUE 
DO  20  J>1#3 
00  20  ialSTART'NTDXM 
TEMP  Cl, J)«OATA(I+NTDIM# J) 

20  CONTINUE 
00  30  jal,3 
00  30  1*1, NTOIM 

PROCCI, J)aTEMP(X, J)*SCPACtOATA(I, J) 

30  CONTINUE 
RETURN 
END 
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SUBROUTINE  RMS 


SU-kOUTI'«t  RMS(PROC,CTIME,TSQD,TlL,TIE,F,TEMP, ISUM) 

Dl'.EMSION  F ( 3) 

01  PENSION  PROC(2l/0,3)»TSQD(i),CTIiiECi),T(3),TEMP(3J 

DO  1 Jsl,3 

F(J)si?,D 

TS'JOC  J)=0.0 

CTIMfc  ( J)=3.1? 

TEmP( jJao.B 
T ( J ) s ft  . 

I F C I S U v . rj t . •)  GO  Tj  5 
NL-1 IL/ a • 

DC  TO  6 

5 NL-TU./2. 

6 NLsnL-2>! 

ANL*NL 

DU  IS  1=1, Nl 
T (7) =T ( J) +P«QC ( l f J) **? 

IS  CUvTlNUE 

IFCT(Jl.GT.U.O)  GO  TO  17 
F( J)=0.2 
TS'JD  ( J)  =0, 0 
CTi«t(j)si.mi 
TEMP( J)3p.0 

*MH I T t C 3#  16) 

16  FORMATCSX, '♦♦♦♦♦  OTA  TRACE  HAS  BEEN  ZEROED  ♦♦♦♦♦') 

GO  TO  1 

17  F(J)s3w<T(T(J)/ANL) 

Ft J)=2?.*alOG10(F(J) ) 

Tr-lP(.JlSJ.^ 

NEsTIt/2. 

IF ( ISUO,  Jt.. ) GO  TO  18 
ML  = 2*'V'L 

16  GO  ?.*  I = ’:L^E 

IF(ABS(PROC(I, J)) .GT.ABS(TEMPCJ)J)  GO  TO  10 
Go  TO  20 
10  10=1 

TEMP(J)=PROC(I,  J) 

23  CONTINUE 

CTIME(J)52*I0 
TS3D(J)=AiiS(TEMP(J)  ) 

TSy'0(J)=2-'.*AL0Gl7i(  rSOO(J))-F(J) 

1 CONTINUE 

return 

Eng 


SUBROUTINE  LOVFIL 


r 

C 

c 


c 

r 

C 


c 

c 

c 


SUBROUTINE  10VFILCPPUF, IPTSEG, ANG8ET) 

DIMENSION  Pbllr  (64,3,2) 

R R S h £ r u . 

TPSMI=4.4 

TlSnls^i.W 
dE T A = n , 3 

Trrot  = ^.,; 

T I *t  0 T s t . iti 
7 RUNShsu,  .0 
T IU'-SHs  '/ 

HTSEGbIPISEG 
CCjNVi^,si8r*,P/3.  U 159 
DU  6«>  1 = 3,  IPTStG,  2 
IF(Pbi.iF(I  + l,l,l))10,2Pf  1« 

14  PMl  = AT»-*2(PBUF(Itl,  1,  1)  ,PBUF(1, 1 , 1) ) 

GO  TD  i ) 

20  ohIs^.v- 

SiIFl  Tint  L-fISiN 

34  RRbPlaPaJP  Cl,  3,1  )«COS  (PHI)  +PBUF(I+  1,3,1  )*S  IN  (PHI) 
r*Shl3»KUF(I,2f  U*COS( PHI) 4>PBUF(  1+1,2,  1)*SXN(PHX) 
TISrtlsR  ii)F(  1 + 1,2,  1)*CCJ>IPH1  J-PBUFC  1,2, 1)*S  IN  (PHI) 
RISHIs33UF(Iti,3, 1)*C0S(PHI)«PBUF(I,3,1)*SIN(PHI) 

COMPUTE  IK  TEST  SETA 

6lT4sA  f A'i(SAKG) 

DEGGt r=3ETA*C0NVHU 
A“36ETsA3S(OEG3ET) 

IF  ( A*SiiET.GT.  ANQRET)  GO  TO  40 

>01  ATE  8Y  MInuS  bETA 

T4PP TsT«SHl*C03( BETA )-RKSHl*S IN (BETA) 
TlHOT=TI:j*I*COS(META)-mShI*SlN(BETA) 

Si-'lKl  HACK  TO  ORIGINAL  TIME  ORIGIN 


I‘5u  '.SHsrRROT*CGS(PHl)-TIROT*SIU(PHI) 
TTU  .bH3rw*0T*SlU(PHn+TIR0T*C0S(PHI) 
GO  TO  53 
4 o T Ru  ‘Jo  Hr  1 # U 
T 1 L N 3 H s . 1 . .1 

5*  Pt"JF(I,2,  lJsTWiJNSH 
P?LF(l+t,2,  lJsTIi  NSH 

60  CONTINUE 
RETURN 


SUBROUTINE  RAYFIL 


SUritfO'Jri\E  kAYFILCP&UF#  IPTSEG,  IP0*ER) 

DIMENSION  PPUF(64,3#2) 

RVPO*=3.? 

R P P 0 ' s 10  • % 

« <f  rCTs/l.^ 

RM  TOT  314,1 

DO  5^  1*1,  IPTSEG, 2 

CPRWs»-3'JF(  I,  i,  1)*P8UFCI,3,  1 ) tPBUF  ( !♦  I , l , 1 ) *PBUF  C 1 + 1 , 3,  1 ) 
CPkIsP3JF(I,1,  1)*PHUF (1  + 1,3#  D-P8UF (I*  1,1,1  )*PBUF (1,3, 1) 

IFCCF-RI.NE.P.tf)  GO  TO  2P 
MFILs'-’.  A 

GO  10  3-1 

2v)  RFIL=CRlI/SimT(CPHI»r*?TCPRR**2) 

3k)  DC  4’/.  J=1,IP0*ER 

IF(RFIL.LT..0C?71)  RFILS0.0 

4.)  RFIl=PFlL**2 

PdUF(I, l, l)sPBUF(I, It 1)*«FIL 
PHUFdfl,  1,  l)*PBUFdM»  1,1)*RFIL 
PttUF  ( I » 3,  l ) sPB'JF  d » 3,  1 ) *RFXL 
PbUFtI+1, 3,  n«PBUF(I*l,3#l)*HFIL 
SO  CONTINUE 
RETURfv 
EnC 


SUBROUTINE  FILTBP 


SimOuTIHG  FILT3P(FIN,MTDIM,FL0,FHI, TAPP) 
DIMENSION  FINC210(!,3J  ,FILT(2108} 

Entoi-s  .toi* 

0ELTF=1  ,?!/(P,i?*ENTi)IM) 

FSTs( (FL0-TAPK)/DELTF)*2,H 
FL0Us(Fl0/jFLTF)*2.i<} 

Fhii  = (Fh1/[;elTF)*2.(? 

KST=Fsr 
KLOOsF L jO 
KMIIsHlI 

TA:mAlF=1.'V(F^OO-FST) 

JPTS=KLOJ-KST 
KLt  '■  sK-tl  I-t\LOO 
On  1-1^  I = 1,nTOIm 
F 1 1_ T (I)s«3.2 
lfel  COW  I ioE 

DO  Hi  I si, JPTS 

A 1 s I 

F I >.  T CI  + KST)sAl*TANALF 
in  COWIvjE 

DO  12.’  i*l,KLEN 
FILT  (•'L’lO+I  )sl,t 
12.1  CON  T I <,'JE 

A JP I 3s JPTS 
DP  13’  tsl.JPTS 
A 1 = I 

FTLT ( KnI 1+ I ) s A JPTS- A I ) * T ANALF 
13.3  continue 

DO  Jr  1,3 

CO  2fU?  lsl,NTDIM 

FInCI,  J)sFIfJ(I,  J1*FILT(I) 

2:30  COnTI*!I)L 
WETUPn 
E NO 
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SUBROUTINE  PLTRTN 


SUBROUTINE  PLTRTN(OATA^NPTS»T,TL»TR/TE»SF) 
OImENSIuN  OATA(MPTS) 

ENsnPTS 

AlFNGT=(EN*T)/130.0 

TlL=(-ALF.NGTfTL*.«?n 

TIn=C-ALENGT+T«*,31) 

TIE=(-ALENGT+TE*.P1) 

KALtriGsALEMGT 

DIST=*ALENG 

FlNLEMsAuENGT-OlST 

BLENGTs-ALENGT 

OXa-*i,il*T 

S = 3.5> 

00  5 Isi#f4pT5 

S IF(AaS(JATA(Z)).GT.S)  SsARS ( DAT  A ( I ) ) 

IF  (S.tO.0.0)  G(.  TO  7 
SF's  !.  VS 
GO  Tu  8 

7 3 F s 

8 CALL  PLOTO.O,-  *.6»-2) 

DO  1 *C<-l»KALEr-t; 

CALL  PLOT  (1. 3, 0.3, 2) 

CALL  PLOT,  (l.fj,  *.  !,  2) 

CALL  PLOT  Cl.>V-'3.1r2) 

1 CALL  plot  ( l.ti,  *.11,-2) 

TALL  PLOT  (FINl.EN,0.t<i,-n 
CALL  plOT(TIE,;<.0,3) 

CALL  PLOTiTK#- O.A,2) 

CALL  PLOTCTIP,.  . Jr  3) 

CALL  PLUT(TIHr“0.6r2) 

CALL  PLOT  ( T iLr  K r 3 ) 

CALL  PLOT (TILr“3.br2) 

CALL  P , 0 T ( T I L r r 2 ) 

E>s".0 

EYE  = DATA(  IPTS)*SF 
CALL  PlDT (DXr  £' Er  3) 

DO  2 -IslrNPTS 

tAstXf.'V 

JJs'iPTS-J+1 

EYsOATA(JJ)*SF 

2 CALL  PLOT  CEXrtY,2) 

CAll  PlOT  (oLEr‘GTr3.3r3) 

CALL  PLOT  (flLE^GT,-0.8r-2) 

RETURN 

EnO 


APPENDIX  B 
DETECTION  STATISTICS 


The  figures  of  this  appendix  give  the  detection  statistics  for 
southern  Kamchatka  (SKAM)  single-site  data,  central  Asia  (CENA)  single- 
site data,  and  central  Asia  (CENA)  beam  data.  The  upper  portion  of  each 
figure  shows  a histogram  indicating  the  number  of  detected  and  non-detected 
events  for  each  m^  increment.  The  lower  portion  of  each  figure  shows  the 
maximum  likelihood  detectability  curve  for  the  statistics  represented  by  the 
histogram. 
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DETECTION  STATISTICS  FOR  SKAM  BANDPASS  FILTERED 
LQ-T  COMPONENT  - SINGLE -SITE  DATA 
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DETECTION  STATISTICS  FOR  SKAM  ORIGINAL  TCA  PROCESSED 
LR-V  COMPONENT  - SINGLE-SITE  DATA 
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FIGURE  B-4 
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DETECTION  STATISTICS  FOR  SKAM  NEW  TCA  PROCESSED 
LQ-T  COMPONENT  - SINGLE-SITE  DATA 
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DETECTION  STATISTICS  FOR  CENA  BANDPASS  FILTERED 
T.R -V  COMPONENT  - STNfTl . F! -RTTF  DATA 
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DETECTION  STATISTICS  FOR  CENA  ORIGINAL  TCA  PROCESSED 
LQ-T  COMPONENT  - SINGLE-SITE  DATA 
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DETECTION  STATISTICS  FOR  CENA  NEW  TCA  PROCESSED 
LQ-T  COMPONENT  - SINGLE-SITE  DATA 
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DETECTION  STATISTICS  FOR  CENA  BANDPASS  FILTERED 
LR-V  COMPONENT  - BEAM  DATA 
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